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THE CENTERED ONE -DIMENSIONAL UNSTEADY EXPANSION 
OF A VIBRATIONALLY RELAXING 
NITROGEN-OXYGEN MIXTURE 

By Laurence N. Connor, Jr., and Frances W. Taylor 
Langley Research Center 

SUMMARY 

The one -dimensional unsteady expansion of a mixture of gases is analyzed for con- 
ditions at which the deexcitation of molecular vibrational energies is the predominant 
rate process. Methods of calculation are formulated for treating the expansion as either 
an equilibrium or a nonequilibrium process. In the nonequilibrium analysis the method 
of characteristics is applied in a Lagrangian frame of reference, and the simultaneous 
vibrational deexcitation of a number of molecules can be handled. 

Calculated results for equilibrium and nonequilibrium expansions are presented for 
a range of initial temperatures from 2000° K to 4000° K. The velocity, thermodynamic- 
property, and vibrational-energy variations are presented along specific particle paths 
designated by the scaling parameter. The general trends seen in the results agree well 
with trends expected from physical reasoning and are very similar to those exhibited by 
vibrationally relaxing flow in a hypersonic nozzle. Curves are presented to permit prac- 
tical applications of the calculated results. 

INTRODUCTION 

The centered one -dimensional unsteady expansion has recently assumed new impor- 
tance as a process for producing high-enthalpy flow for rate measurement experiments 
and reentry simulation. The expansion tube proposed in reference 1 is a prime example 
of the potential use of such an expansion to generate high-velocity flow without subjecting 
the gas to extreme stagnation temperatures and pressures. 

The rapid expansion from elevated temperatures experienced by a gas mixture in an 
unsteady expansion necessitates an analysis which includes the effects of the nonequilib- 
rium processes which exist in the flow. The present investigation is intended to present 
such an analysis for a limited range of conditions. For the temperatures and pressures 
considered, the nonequilibrium processes considered can be restricted to those governing 
the relaxation of the vibrational energies of the molecular species in the gas mixture. A 


constant-composition gas mixture in which the vibrational energies are rate controlled 
is used in the analysis. 

The structure of a centered unsteady expansion in a vibrationally relaxing gas is 
first analyzed in reference 2. In reference 3 a more realistic vibrational rate equation 
is used to treat the same problem. 

The purpose of the present work is to extend the existing knowledge of vibrationally 
relaxing unsteady expansions by treating a mixture of gases and the associated rate equa- 
tions. The method of characteristics is used and calculations are performed in a 
Lagrangian coordinate system. Results are presented for a nitrogen-oxygen mixture 
over a range of conditions in which vibrational relaxation is the principal nonequilibrium 
process. The results given should prove valuable in the design and analysis of facilities 
operating in the range considered. 


SYMBOLS 

Ai,Bj,Ci constants used in vibrational-relaxation-time expressions 

a speed of sound 

b Lagrangian coordinate denoting a fixed mass element 

Cp heat capacity at constant pressure 

f^ factor used to provide exclusion and inclusion of rotational energy of 

monatomic and diatomic molecules, respectively (fi = 0 for mon- 
atomic molecule; fi = 1 for diatomic molecule) 

h enthalpy per unit mass 

i species identification 

l distance from center of expansion 

n total number of species 

p pressure 

R universal gas constant 
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s 

T 

t 

U 

u 

X 


y 

y 

M 


H 


P 


Subscripts: 

A,B,D,P 

AP,BP,DP 

c 


specific entropy 

temperature 

time 


piston velocity 


velocity 


position coordinate of fluid particle 
mass fraction of ith species 
Lagrangian coordinate divided by t 
specific heat ratio 

characteristic vibrational temperature of ith species 

molecular weight 

molecular weight of ith species 

density 

vibrational energy of ith species and gq^ represent the vibra- 

tional energies of N 2 and O 2 , respectively^ 

vibrational relaxation time of ith 




points in figure 2 

denote average of values at two points in figure 2 
required for complete expansion 
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equilibrium conditions 
results for frozen flow 
evaluated at intermediate condition 
initial conditions of undisturbed gas 
evaluated with pressure held constant 
evaluated with specific entropy held constant 
evaluated with density held constant 

evaluated with vibrational energy of ith species held constant 
denote the three characteristic directions 

A bar over a symbol denotes a nondimensional quantity. 

THEORY 

Physical Model of the Expansion 

A physical picture is necessary in order to understand the unsteady one-dimensional 
expansion and its related boundary conditions. The model used in the present work is 
illustrated in figure 1. 

A semi-infinite reservoir of constant cross-sectional area containing a gas mixture 
in a state of thermodynamic equilibrium is considered. At time zero, the piston is sud- 
denly withdrawn with a constant velocity. A rarefaction wave which propagates into the 
undisturbed gas in the reservoir is produced. The first infinitesimal pressure pulse 
propagates into the undisturbed gas with the velocity of the speed of sound of that gas 
(-a Q ). This pulse is followed by an infinite number of such pulses traveling at the speed 
of sound of the gas into which they are progressing (u n - an). These infinitesimal pulses 
continue until the properties at a point in question are identical with those at the piston 
face. Figure 1 graphically depicts the structure of the expansion in distance-time coor- 
dinates. For a perfect gas, the centered rays of the expansion are straight lines along 
which the flow properties are constant. Typical particle paths are indicated on the x - t 
diagram by dashed lines. 
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The temperature and pressure of the undisturbed gas are assumed to be high and 
the vibrational energies of the molecular species of that gas are initially in an excited 
state. The temperature and pressure decrease as the molecules are processed by the 
expansion, and vibrational deexcitation must occur if thermodynamic equilibrium is to be 
maintained. However, the changes in vibrational energies are rate controlled and there- 
fore rate processes must be included in the analysis of the expansion. 

An indication of the effects of rate processes can be seen by considering the various 
particle paths shown in figure 1. Path I represents a particle which was very near the 
piston face initially. This particle has very little residence time in the expansion fan and 
thus the vibrational energies have insufficient time for adjustment to equilibrium. If very 
little adjustment occurs, the particle essentially experiences a frozen expansion and the 
vibrational energy is not converted to the translational and rotational energy modes. 
Particle in, on the other hand, has a relatively long residence time in the expansion and 
the rate processes are permitted adequate time for the flow to relax to equilibrium. 

Such a particle undergoes an equilibrium expansion. The energy contained in the excited 
vibrational states returns to the translational and rotational modes. Between the two 
limiting cases of frozen and equilibrium flow is the nonequilibrium case, illustrated by 
particle path II. The residence time is of a magnitude that allows some relaxation of the 
rate -controlled vibrational energies toward equilibrium. Only a portion of the vibrational 
energies is returned to the flow. The principal portion of this investigation concerns the 
nonequilibrium region. Solutions for equilibrium and frozen flows are also presented 
for comparison purposes. 


Governing Equations 

Flow equations .- The analysis is formulated in a Lagrangian frame of reference to 
facilitate the establishment of particle paths along which the rate equation can be applied. 
In Lagrangian coordinates, attention is given to what happens to an individual fluid parti- 
cle in the course of time. Each individual fluid particle must therefore be labeled. In 
the present work, this labeling is accomplished by designating each particle according to 
its position coordinate x in the x - t plane at some reference time t = t 0 . This value 
is called b, the Lagrangian coordinate specifying the particle being considered. Inde- 
pendent variables in the Lagrangian coordinate system are therefore b and t. The 
governing equations for a one -dimensional flow in such a system are presented in refer- 
ence 4. These equations, expressing the conservation of mass, momentum, and energy 
in a Lagrangian frame of reference, can be written as 


+ o 

ab p2 at 


( 1 ) 
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( 2 ) - 


and 


, ®3i + ^ = 0 

0 9t 8b 


9h 1 9p _ q 
9t " P 9t 


(3) 


provided body forces, viscous forces, diffusion, and heat conduction are neglected. 

State equations .- If dissociation, ionization, and electronic excitation are neglected, 
the caloric equation of state can be expressed in the form 


h = 



(4) 


The exclusion and inclusion of rotational energy of monatomic and diatomic molecules, 
respectively, is provided by use of fi. For odd values of i, fi = 0; for even values 
of i, fi = 1. In the present application, the mass fractions of the various constituents 
of the gas mixture are assumed to remain constant throughout the expansion. 

For the nonequilibrium analysis, the vibrational energies are rate controlled and 

are thus dependent upon the process which the gas has undergone. Equilibrium vibra- 
tional energies are calculated from the expression 



(i = 2,4,. . .,n) (5) 


which can be derived from quantum mechanical considerations by using the simple har- 
monic oscillator model for diatomic molecules. 

The gas is assumed to be a mixture of perfect gases. It therefore obeys the 
thermal equation of state 


where 




and is constant for a given gas composition. 


( 6 ) 
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Rate equations.- When the residence time in the expansion is not adequate for the 
vibrational energy of a molecule to relax to its equilibrium value, a rate equation must 
be used to describe the relaxation process. A form of the vibrational rate equation 
which indicates that a whole assembly of simple harmonic oscillators in various vibra- 
tional energy levels acts in such a way that a single relaxation time exists for the total 
vibrational energy was first established in reference 5. Such a form of the vibrational 
rate equation written in Lagrangian coordinates is 


30^ _ °i,eq " °i 
0t Ti 


The vibrational relaxation time is taken to be of the form 


T i = 


AiT B i exp^CjT" 1 / 3 ) 
— (-?) 


1 - exp 


(7) 


( 8 ) 


The expression describes quite accurately experimentally measured values of vibrational 
relaxation times and is essentially the form predicted by the theoretical work in refer- 
ence 6. The constants used for the various gases in this investigation are discussed 
subsequently when actual calculations are considered. 

Nonequilibrium Solution 

Characteristic equations .- The governing equations of a nonequilibrium one- 
dimensional unsteady expansion - equations (1), (2), (3), and (7) - constitute a system 
of quasi-linear partial differential equations which is hyperbolic in nature. Since real 
valued characteristics exist for such a system, the method of characteristics is used for 
obtaining a solution. The characteristics are lines across which the derivatives of the 
dependent variables may be discontinuous, and the slopes are found by equating the 
matrix of coefficients of the set of governing equations to zero and evaluating the 
resulting expression. 

In the Lagrangian frame of reference the + and - characteristic directions are 
found to be 



(9) 
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Along these characteristic lines, the respective compatibility relations 



exist. It should be noted that a^ refers to the local frozen speed of sound which is 
defined as 


H = 



( 11 ) 


for the gas model used herein. 

The third characteristic direction is identical with a particle path and is written as 



Along this line, both the energy equation 

dh 1 dp _ q 
dt ' P dt 

and the rate equation 

dq gj,eq - °i 

dt ** Ti 


( 12 ) 


(13) 


(14) 


must be satisfied. 

Modified Lagrangian coordinates .- In order to simplify the application of the bound- 
ary conditions to the problem, one further modification of the independent coordinates is 
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necessary. This modification is made by dividing the Lagrangian particle coordinate b 
by the time variable t to define a new variable y. The independent variables then 
become y and t, where y is defined by the relation 

y = * (15) 

The transformation maps the corner of the expansion fan into a straight line and readily 
permits application of the boundary conditions for t — 0. In addition, the grid network 
becomes more nearly rectangular in the y - t plane. 

In the modified coordinates, the equations which prescribe the characteristic 
curves, equations (9) and (12), become 


and 



(16) 


(17) 


respectively. The compatibility relations remain unchanged. 

Nondimensionalization of the characteristic equations .- The following dimension- 
less quantities are defined: 


3- 

H = ~ — 

p = £- 

h-A 1 

^,0 

Po 

h o 

n-JL 

a - &e< l 

U - JL 

Po 

a eq ~ a 

a eq,o 

Uc 

b = A 

- ^ 

y = JL_ 

bo . 

b 0 

H, 0 


(18) 


The subscript o refers to the value of the quantity in the undisturbed gas before the 
piston is withdrawn. The velocity U c is that required for a complete frozen expansion 

to p = 0. It can be readily shown that U c = — ^2 The reference length b 0 is 

the initial x position of some reference particle. 
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The substitution of these expressions into the previously developed equations 
yields the dimensionless characteristic equations for the modified Lagrangian coordi- 
nate system. Along the + and - characteristic curves 



j 


the ordinary differential equations 
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dp 




(19) 


( 20 ) 


must be satisfied; along the particle path 



( 21 ) 


the energy and rate expressions become 


dh_Jl2_ld£ = o (22) 

dt h 0 p Q p dt 

and 

do i _ _^o_ q i,eq ~ q i (23) 

dt ^,0 T i 


y — 1 

Scaling law.- For a fixed initial composition, the factor — is constant and 

2y 0 

1/U C and P 0 ^hop 0 are functions of T 0 . Equations (20) and (23) reveal that, for 
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a fixed value of the product p 0 b 0 and a given value of T 0 , the rate equation remains 
the same. Thus, for fixed values of T 0 and the scaling parameter p 0 b Q , the set of 
characteristic equations is invariant and generates identical flow fields. Since for a 
given T 0 , t 0 (the time the reference particle enters the expansion fan) is directly pro- 
portional to b Q , the product p Q t 0 can also be used as the scaling parameter. 

Evaluation of enthalpy derivatives .- In order to apply the characteristic equations, 

it is necessary to obtain expressions for , ( jp] , and aj in terms of p, p, 

V a i/p,p \ P 'p,aj 

and Oj. By utilizing equations (4) and (6) and remembering that x^ is constant, it 
becomes apparent that 



Similarly, the frozen sound speed can be shown to be 



(24) 


(25) 


(26) 


Boundary conditions .- Two basic boundary conditions are prescribed for the expan- 
sion region. First, the undisturbed gas in the reservoir is assumed to be in a state of 
thermodynamic equilibrium - that is, at t = 0 before the piston is withdrawn, the equi- 
librium properties of the gas mixture are known. This condition may be written for the 
y,t coordinate system as 


P(M) = P r 


P(l,t) = P f 


u(l,t) = u 0 ^(1,0 = <T i>0 


(27) 


11 


The second boundary condition is obtained by considering the unsteady expansion as t — 0. 
A particle processed by the expansion for this second condition has a residence time in 
the expansion which approaches zero. The particle adjacent to the piston face experiences 
such an expansion but, because of the short residence time, no relaxation of the rate- 
controlled vibrational energies toward equilibrium occurs. It can be verified mathemat- 
ically that, as t — 0, the properties throughout the expansion are those of a frozen expan- 
sion and the vibrational energies remain at their initial values. This boundary condition 
can be written as 


P(y,0) = p f (y) = y y o +1 

2 

p(y,0) = P f ( y) = y y ° +1 > 

7o-l 

u(y,0) = uj(y) = 1 - y 7 ° +1 + u Q 
a i(y»°) = CT i,o 


(28) 


These two boundary conditions are sufficient to allow calculation of the expansion region 
by the method of characteristics. 


Solution Procedure 

Finite-difference equations .- In order to establish a method for obtaining a numeri- 
cal solution, the characteristic and compatibility equations must be written in finite- 
difference form. By using the general grid nomenclature established in figure 2 and the 
finite -difference approximations, equations (19) and (20) may be written as equations (29) 
to (32). Along 
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applies. Along 


yptp - y B t B = 



the compatibility equation 


(31) 


holds. 



The particle path is described by 


(32) 


yptp = y D t D 


(33) 


and, along this path, the energy and rate equations 



must be satisfied. These finite-difference equations can be used to prescribe a solution 
procedure. 

Calculation outline .- The grid network generated by the flow field calculation is 
represented in figure 3. The equilibrium boundary condition prescribes properties along 
the line y = 1.0, and along the line t = 0 the frozen boundary condition applies. There- 
fore, the calculation begins in the lower right-hand corner of the grid network. A 
marching process is used to carry the calculation point by point through the desired flow 
field. 


The key operation in computing the flow field is the calculation of values at any 
unknown point by using the known values at three adjacent points and the finite-difference 
characteristic equations. This operation can be performed by an iterative process, the 
steps of which are outlined herein. The nomenclature used is that defined in figure 2. 
The objective is the determination of the coordinates, properties, vibrational energies. 
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and velocity at point P. Those quantities are known at points A, B, and C from the bound- 
ary conditions or preceding calculations. 

As the first assumption, the values of the properties at point P are assumed to be 
identical with those at point A. For equilibrium flow the expansion is a simple wave 
region and properties along a characteristic line similar to line AP are constant. 

Changes in properties along line AP in the present study are due only to nonequilibrium 
effects and are generally smaller than those along line BP which are basically due to the 
expansion process. Equating property values at point P to those at point A should there- 
fore provide more rapid iteration convergence than equating them to the average of values 
at points A and B. Equations (29) and (31) are solved simultaneously to give the coordi- 
nates of point P. The compatibility relations along lines AP and BP, equations (30) and 
(32), are applied to obtain improved values of velocity and pressure at point P. By using 
equation (33) the particle path through point P is extended backward and its point of inter- 
section, point D, with line CA or CB is obtained. Properties at point D are determined by 
using a linear interpolation between points C and A or C and B. Vibrational rate equa- 
tions of the form given in equation (35) are applied along the particle path (line DP) to 
obtain the vibrational energy value at point P. The energy equation, equation (34), is also 
used along the particle path to get an improved value of enthalpy. The caloric and thermal 
equations of state yield new values of temperature and density at point P. Thus, an 
entirely new set of quantities is calculated for the point. These improved values at 
point P are then used in a second-order approximation scheme for the next iteration. 

When two consecutive iterations produce temperature values within an acceptable accu- 
racy, the calculation is performed at the next point. The process is repeated until values 
have been calculated at all desired points. 

Computer program .- A program, written in FORTRAN IV language for the IBM 7094 
electronic data processing system, was developed to perform the required calculations. 
This program may be used to treat a number of chemical species and their associated 
vibrational rate equations. Numerical accuracy of the program is good provided that the 
pressure drop throughout the expansion does not exceed two orders of magnitude. In the 
present application, this limitation is not serious inasmuch as the vibrational energies 
have either essentially frozen or become insignificant by the time the gas has experienced 
that degree of expansion. 

Calculated results may be printed out in two ways. The properties, velocities, and 
vibrational energies can be printed out either at each point in the flow field for the coordi- 
nates of that point or at various points along specified particle paths. A numerical inte- 
gration routine is included in the program to calculate the distance traveled by the specific 
particles. 
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* Equilibrium Solution 

When the residence time of a gas particle in the expansion is adequate to allow 
adjustment of the vibrational energies to their equilibrium values, the expansion is isen- 
tropic. The analysis of an equilibrium expansion is presented herein. The composition 
of the gas mixture is assumed to remain constant throughout the expansion. 

Since the process is isentropic, the expression 

ds = %f!I-£&=0 (36) 

T |i P ' 

applies and can be integrated to give 



The differential equation describing the flow through an isentropic unsteady expansion 


du = 



/ r» n\ 
\ 00 ) 


is numerically integrated to obtain the velocity throughout the expansion. The equilibrium 
sound speed a eq can be shown to be 
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for the gas model in the present study. 

Values of T/T 0 are chosen in small increments from 1.0 to zero. For each value 
of T/To, the pressure is calculated by using equation (37) and the density is found by 
using the thermal equation of state. The value of the coordinate y is found for each 
value of T/T 0 from the expression 


y = P^eq 


(40) 


Equation (38) is numerically integrated over the expansion to obtain values of velocity. 

Description of Gas Model 

A gas mixture of 21 percent O 2 and 79 percent N 2 was used to approximate air for 
calculation purposes. Conditions at which vibrational deexcitation is the dominant rate 
process were chosen. Dissociation of molecular species is assumed to be negligible. 
This assumption is made valid by the prescribing of minimum pressures for each initial 
temperature. The presence of nitric oxide NO is ignored, as it becomes appreciable 
only at the higher temperatures considered and, at those temperatures, produces very 
small changes in the thermodynamic properties. 

The constants used in computing the vibrational relaxation times of O 2 and N 2 were 
obtained by curve fitting the values of reference 7 and are listed in the following table: 
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Physical constant 

Value for - 

02 

n 2 

0i, °K 

2270 

3390 

Ai, atm -sec 

0.1424 

560.61 

Bi 

0.3718 

-0.0995 

Ci 

179.33 

163.63 


Experiments described in reference 8 indicate that N 2 vibrational relaxation times in 
expanding flows are approximately 15 times as fast as those measured in compressive- 
flow experiments such as that of reference 7. These findings are substantiated in refer- 
ence 9 for the expansion of both pure N 2 and air. The present calculations were made by 
using data from reference 7 as a base. By changing the scaling parameter p 0 t 0 for a 
particular expansion, the results can readily be adjusted so that the faster deexcitation 
rates of references 8 and 9 may be used. The procedure for making this change is 
described subsequently. 

It is further assumed that the N 2 and O 2 vibrational relaxation processes are inde- 
pendent and therefore no exchange of vibrational energy between the separate species 
occurs during collisions. The effect of an exchange in a compressive-flow situation is 
discussed in reference 10 where it was concluded that the vibrational relaxation rate of 
the N 2 molecule is increased by a vibrational energy exchange with the more rapidly 
relaxing O 2 molecule. No experimental evidence of a similar exchange in expanding 
flows is available, and vibrational energy exchange between excited molecules of the 
separate species is therefore neglected in the present analysis. 

RESULTS AND DISCUSSION 

The structure of a centered one-dimensional unsteady expansion in an N 2 and O 2 
mixture closely approximating air was calculated for a number of flows. Results are 
presented along particle paths through the expansion. The correlating parameter p 0 t Q , 
the product of initial pressure and the time the particle entered the expansion, designates 
the various particle paths. Numerical results are presented in figures 4 to 7 for initial 
temperatures from 2000° K to 4000° K in 500° increments. The pressure limitations 
within which the results are valid are stated in the figures. These limits are defined as 
the pressures at which the vibrational energy is one order of magnitude larger than the 
energies involved in dissociation. Vibrational energies, thermodynamic properties, and 
velocities are displayed as functions of the reduced time coordinate t/t 0 , where t is 
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the actual particle time and t 0 is the time at which the particle entered the expansion 
fan. 

As already mentioned, recent experiments have indicated that vibrational relaxation 
rates in expanding N2 and air are about 15 times as fast as the values measured behind 
shock waves. In the present calculations, the shock data of reference 7 are used as a 
base. In order to adjust the calculations to allow for the assumption of the more rapid 
rates, the scaling parameter p 0 t 0 appearing in figures 4 to 9 must be multiplied by 
1/15. 

The trends exhibited in the unsteady expansion of a vibrationally relaxing gas are 
found to be quite similar to those found in the quasi-one-dimensional nonequilibrium 
expansion in a hypersonic nozzle. These similarities are pointed out in the ensuing 
discussion. 

Figure 4 depicts the variation of the vibrational energies of particles processed by 
an unsteady expansion. The equilibrium and frozen curves represent two limiting cases. 
The intermediate curves indicate the amount of vibrational nonequilibrium which exists 
along various particle paths. 

Several observations can be made. First, the total amount of vibrational energy 
converted back to translational and rotational energy modes is highly dependent upon the 
pressure and the residence time of the particle, as characterized by the product p 0 t Q . 

An increase in initial pressure or residence time results in increased adjustment of the 
vibrational energies toward their equilibrium values. 

In some near-equilibrium flows, the O 2 vibrational energy is less than the value 
which would be obtained in an equilibrium expansion. This phenomenon can be seen in 
figure 4(b). The vibrational energy of the O 2 molecules approaches equilibrium much 
more rapidly than the vibrational energy of the N 2 molecules. The O 2 molecules attain 
values of vibrational energy which are in equilibrium with the local translational temper- 
ature. However, since the N 2 vibrational energy is still in a nonequilibrium state higher 
than the equilibrium value, the translational temperature remains below that of an equi- 
librium expansion. The vibrational energy of the O 2 molecules adjusts to the equilibrium 
value for this temperature and therefore undershoots the value which would be attained in 
an expansion where both the N 2 and the O 2 vibrational energies remain in equilibrium 
throughout. 

Another observation of significance is that the vibrational energies relax rapidly 
during the initial portion of the expansion and then the rate of adjustment quickly dimin- 
ishes and the vibrational energies approach constant values as further expansion occurs. 
This freezing phenomenon can be attributed to the fact that the relaxation times are 
strongly dependent upon both temperature and pressure, as is evident from equation (8). 
As the gas is processed by the expansion, the temperature and pressure decrease rapidly 
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and the relaxation times thus increase until the change of vibrational energy becomes 
negligible. Calculations made by carrying the expansions to pressures lower than those 
presented indicate that any further vibrational relaxation is negligible. This freezing of 
vibrational energy along a particle path in an unsteady expansion is very similar to the 
freezing which occurs in hypersonic nozzles. A comparison can readily be made between 
the results presented herein and the results presented in reference 11 for the hypersonic 
nozzle. 

Thermodynamic properties of the expansion also reflect the degree of nonequilib- 
rium which exists in the flow. Static temperature, defined for this work as the transla- 
tional and rotational temperature of the gas, proves to be the property most affected by 
the relaxation process. Vibrational relaxation toward equilibrium returns energy to the 
translational and rotational modes. This energy manifests itself as an increase in the 
temperature and velocity of the gas. If vibrational equilibrium exists, the maximum tem- 
perature exists throughout the expansion. As p Q t 0 increases, the energy returned to 
the translational and rotational modes from the excited vibrational state increases, and 
so the temperatures increase. This behavior is demonstrated in figure 5. The depend- 
ence of pressure and density upon the degree of nonequilibrium is indicated in figure 6. 

The effect of the rate processes on the flow velocity is shown in figure 7. As the 
parameter p 0 t 0 is increased, vibrational relaxation increases, more energy appears in 
the translational mode, and higher velocities are attained. Maximum flow velocity exists 
for equilibrium flow. 

Figures 8 and 9 are presented to permit practical applications of the results. If 
T 0 and the product p 0 t 0 are known, the effects of nonequilibrium on the thermodynamic 
properties and on flow variables can be obtained by using figures 3 to 7. When a desired 
velocity u is reached in the expansion at a distance l from the center of the expansion 
and when the undisturbed equilibrium temperature and pressure are known, figures 8 
and 9 may be used to obtain the value of t 0 . The procedure for obtaining this value is 
outlined in the next paragraph. 

The values of p 0 , T 0 , and l are known, and the reference velocity U c is given 

L . DT 

\ I /'o AV A i 

by U c = 2\l . A value of slightly less than the desired u is assumed, and the 

y p y - 1 

curve for frozen flow in figure 8 is used to get 1/ U c t Q and thus an initial approximation 
for t Q . If this value of t 0 is used to obtain the parameters p 0 t 0 and l /U c t 0 and if 
T 0 is known, the ratio of nonequilibrium velocity to frozen velocity u/uf can be deter- 
mined from figure 9. Since the velocity value used in figure 8 is actually Uf, u is now 
determined. If this u is not the desired velocity, a new value of t Q is assumed and 
the procedure is repeated. This iterative scheme is continued until the desired velocity 
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is obtained. The corresponding t 0 is then the correct one for the nonequilibrium 
expansion and is used to form the parameter p Q t 0 for use in evaluating a given expan- 
sion from figures 3 to 7. 

Values of l / Uct Q greater than those presented in figure 9 represent expansion 
beyond the point at which the flow becomes essentially frozen. For such an expansion, 
accurate approximations of t G can be obtained by using the curve in figure 9 for 

— — =40 in determining u/uf . Also, curves are presented in figure 8 for determining 

u c t G 

t 0 for equilibrium expansions from various initial temperatures. 

CONCLUDING REMARKS 

Methods are presented for the calculation of the one-dimensional unsteady expan- 
sion of a mixture of gases in which the vibrational energies may be either in equilibrium 
or rate controlled. The method of characteristics, formulated in a Lagrangian frame of 
reference, is used in the nonequilibrium analysis. Calculated results, obtained by using 
the IBM 7094 electronic data processing system, are presented over a range of initial 
temperatures from 2000° K to 4000° K. These results, given in terms of flow properties 
and vibrational energies along specific particle paths, agree well with the trends expected 
from physical reasoning. Curves are included to permit practical applications of the 
calculated results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 2, 1966, 

129-01-08-16-23. 


20 


REFERENCES 


1. Trimpi, Robert L.: A Preliminary Theoretical Study of the Expansion Tube, A New 

Device for Producing High-Enthalpy Short-Duration Hypersonic Gas Flows. NASA 
TR R-133, 1962. 

2. Wood, William W.; and Parker, F. R.: Structure of a Centered Rarefaction Wave in 

a Relaxing Gas. Phys. Fluids, vol. 1, no. 3, May- June 1958, pp. 230-241. 

3. Appleton, J. P.: Calculations of the Structure of Unsteady Rarefaction Waves in 

Oxygen/Argon Mixtures, Allowing for Vibrational Relaxation. Tech. Note No. 
Aero.2873, Brit. R.A.E., Feb. 1963. 

4. Courant, R.; and Friedrichs, K. O.: Supersonic Flow and ShockWaves. Interscience 

Publ., Inc. (New York), 1948. 

5. Bethe, H. A.; and Teller, E.: Deviations From Thermal Equilibrium in ShockWaves. 

Rep. No. X-117, Ballistic Res. Lab., Aberdeen Proving Ground, 1945. 

6. Landau, L.; and Teller, E.: Theory of Sound Dispersion. Phys. Zeitschr. der 

Sowjetunion, vol. 10, no. 1, 1936, pp. 34-43. 

7. Blackman, Vernon: Vibrational Relaxation in Oxygen and Nitrogen. J. Fluid Mech., 

vol. 1, pt. 1, May 1956, pp. 61-85. 

8. Hurle, I. R.; Russo, A. L.; and Hall, J. Gordon: Spectroscopic Studies of Vibrational 

Nonequilibrium in Supersonic Nozzle Flows. J. Chem. Phys., vol. 40, no. 8, Apr. 
1964, pp. 2076-2089. 

9. Sebacher, Daniel I.: An Electron Beam Study of Vibrational and Rotational Relaxing 

Flows of Nitrogen and Air. Paper presented at the 1966 Heat Transfer and Fluid 
Mechanics Institute (Santa Clara, Calif.), June 22-24, 1966. 

10. White, Donald R.; and Millikan, Roger C.: Vibrational Relaxation in Air. AIAA J., 

vol. 2, no. 10, Oct. 1964, pp. 1844-1846. 

11. Erickson, Wayne D.: Vibrational -Nonequilibrium Flow of Nitrogen in Hypersonic 

Nozzles. NASA TN D-1810, 1963. 


21 







Dimensionless time 



0 1.0 

Dimensionless modified Lagrangian coordinate, y 

Figure 3.- Complete grid network. 


24 


i 


Nitrogen vibrational energy, <7 , joules per mole Oxygen vibrational energy, a , joules per mole 
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(a) T 0 = 2000° K ; p Q > 0.015 atm. 

Figure 4.- Variation of oxygen and nitrogen vibrational energies in a nonequilibrium unsteady expansion. 


25 






Nitrogen vibrational energy, <7 , joules per mole 0x ^en vibrational energy, or , joules per mole 


| 

! 

i 



Reduced time coordinate , -p 


(M 


0 


10 20 30 40 50 60 70 80 90 100 

Reduced time coordinate, 

l o 



(c) T 0 = 3000° K ; p 0 > 90 atm. 
Figure 4.- Continued. 
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(d) T 0 = 3500° K ; p 0 > 1000 atm. 
Figure 4.- Continued. 








(a) T 0 = 2000° K; p 0 > 0.015 atm. 

Figure 5.- Temperature variation in a nonequitibrium unsteady expansion. 








32 


Reduced time coordinate 
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Reduced time coordinate 
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Reduced time ccx>rdlnate 






Reduced time coordinate, — 
t o 

(0 T 0 = 3000° K; p 0 > 90 atm. 
Figure 6.- Continued. 
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(a) T 0 = 2000° K ; p 0 > 0.015 atm. 

Figure 7.- Velocity variation in a nonequilibrium unsteady expansion. 
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Figure 7.- Concluded. 




1.12 



46 


NASA- Langley, 1967 


L-4788 


Figure 9.- Nonequilibrium-velocity correction as a function of Tq, p 0 to and 
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